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Abstract. We propose a simple mechanical model for relativistic explosions with both forward and 
reverse shocks, which allows one to do fast calculations of GRB afterglow. The blast wave in the 
model is governed by pressures Pp and Pp at the forward and reverse shocks. We show that the 
simplest assumption Pp = Pp is in general inconsistent with energy conservation law. The model 
is applied to GRBs with non-uniform ejecta. Such "stratified fireballs" are likely to emerge with a 
monotonic velocity profile after an internal-shock stage. We calculate the early afterglow emission 
expected from stratified fireballs. 



INTRODUCTION 



Relativistic blast waves from GRBs are believed to produce the observed afterglow 
emission of the bursts. The blast-wave structure is schematically shown in Figure 1. 
To a reasonably good approximation, the whole blast [the region between FS and RS, 
including the contact discontinuity (CD)] may be assumed to have a common Lorentz 
factor Ys \M- The Lorentz factors of the forward shock (FS) and reverse shock (RS) are 
denoted by Yf an d Yr> respectively. Pressures Pp and Pr behind the forward and reverse 
shocks are calculated from the jump conditions. 
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FIGURE 1. An illustrative diagram of 4 regions in a blast wave. Here n\ and «4 are the number densities 
of the ambient medium and the ejecta, respectively. 



STRATIFIED FIREBALL 



We focus here on the early stages of the GRB explosion when both forward and reverse 
shocks are expected to produce significant emission. Our calculations assume spherical 
symmetry. However, the results also apply to beamed explosions at early stages, when 
Ys > Qjet where 0j et is the opening angle of the ejecta. 

The reverse shock (RS) is expected to propagate in the ejecta and produce significant 
emission at radii r ~ 10 16 — 10 17 cm. The RS emission depends on the structure of the 
ejecta that forms at smaller radii r < 10 16 cm, when internal shocks take place Jj]]. 
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FIGURE 2. Lorentz factor (Left) and kinetic luminosity (Right) of ejecta 

Theoretically, one expects the ejecta to emerge with a monotonic velocity profile 
(internal shocks tend to make the profile monotonic). Therefore we study here radially 
stratified ejecta and calculate the emission expected from such "stratified fireballs." 

A simple example model considered here assumes ejecta with radial profile of Lorentz 
factor T e j(x) = (fi - Ji)e~ a% l %b + y 2 shown in Figure 2 (j\ = 500, y 2 = 100, a = 3, 
and Tfc = 30 s is the burst duration). Here x is the Lagrangian coordinate that labels 
the shells in the ejecta by their ejection times. The luminosity of the ejecta L(t) = 
Anr 2 v e j T 2 e j ri4 m p c 2 is also shown in Figure 2. Specifically, L(t) =Lq/(1+ exp[(r — 

Tfc)/T^]) is assumed here with Lq = 10 52 erg/s and % d = 1 s for the timescale of decay. 
The total energy of the burst is given by E e \ = LqTj ln( 1 + e Th ^ Td ) . 



MODELING THE BLAST WAVE 



1. Customary approximation: P F = Pi 
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The customary approximate model blast wave assumes equal pressures at the forward 
and reverse shocks: Pp = Pr JH. This assumption, as we show in an accompanying 
paper, enables a simple analytical solution for y s and pressure P = Pp = Pr: 
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P = -{tf-l)n l m p c 2 



It turns out however that the assumption of equal pressure leads to a contradiction 
with the law of energy conservation. The total energy of the blast wave and ejecta at a 
time t may be found by integrating the 00 component of the stress-energy tensor over 
volume, E = J ToodV. The result is shown in Figure 3. E decreases with time (dashed 
curve in the left panel), i.e. the energy is not conserved. We therefore propose another 
dynamical description. 



2. Consistent mechanical model 

Instead of finding y s from Pp = Pr we solve the differential equation for y s : 

where H = J(U +P) dr is the integrated enthalpy of the blast. The first term on the right 
describes the external force applied to the blast due to the difference between Pp and 
Pr, and the second term describes the adiabatic acceleration of the blast. This equation 
is derived from the exact relativistic fluid equation. It is obtained by integrating over the 
blast and may be viewed as a mechanical equation. The mechanical model is consistent 
with the energy and momentum conservation (solid curve in the left panel of Fig. 3). 
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FIGURE 3. Left: Dashed curve shows the total energy of the blast wave calculated using Pf = Pr 
assumption, solid curve - using our mechanical model. Right: Lorentz factor 743 of region 4 relative to 
region 3 for the stratified (solid curve) and non-stratified (dashed curve) fireballs; our mechanical model 
is used in the calculations. 



3. Light curves from stratified fireballs 

We calculate the afterglow emission using the standard simple model of the postshock 
plasma. It assumes that a fraction e e ~ 0. 1 of the shock energy goes to the electrons 
and accelerates them with a power-law energy spectrum dN/dE e <x EJ P . The postshock 
magnetic field is parameterized by the "equipartition" parameter £g = (B 2 /S7t)U th l < 1 
where U± is the thermal energy density of the shocked medium. 



We discretize the external medium and the ejecta into spherical mass shells m, and 
use the Lagrangian description of the blast wave. Each m,- is impulsively heated at some 
point by a shock front (forward or reverse). We track the subsequent evolution of the 
shell and calculate its synchrotron emission in the shell comoving frame. 

The shock-heated shell suffers radiative and adiabatic energy losses in the expanding 
blast wave, and its magnetic field decreases, so the shell emission weakens with time. 
We find its contribution to observed radiation at a given observed time t t, s . Integration 
over all m ; - gives the exact total spectrum received from the blast wave at a given t \, s . 




FIGURE 4. Left: Optical light curve (R-band) from stratified (solid) and non-stratified (dashed) 
fireballs. The explosion parameters are Eb — 10~ 3 , £ e = 0.1, p = 2.5, n\ = 10 cm~ 3 , and redshift z = 1. 
Right: IR light curve (v b s = 10 13 Hz) from the same models. 

The optical and infrared light curves for two example models (stratified and non- 
stratified) are shown in Figure 4. We find weak optical emission from the reverse shock. 
It is weak because the RS Lorentz factor relative to the blast, 743, is modest and the 
emitted spectrum peaks in IR rather than in optical band. The peak frequency is 

v m ^10 6 fl^Hz, rm= i + ^^e e ( r43 _i). 

p — lm e 

The evolution of 743 is shown in Figure 3 (right panel); its typical value is 1.2 and the 
factor (743 — l) 2 that enters V m is small, ~ 1 /20, shifting the peak frequency to IR. 

The behavior of the RS Lorentz factor 743 in the stratified fireball is significantly 
different from the non- stratified case: it reached maximum and then minimum (Fig. 3). 
Therefore, the RS emission in the stratified fireball is weaker and its peak in the light 
curve is less pronounced, with a smoother rise and slower decay. 
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